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The results of computer simulations of light scattering by a random array of small particles is 
presented. A tensor Green's function method is employed. Results are given for arrays of particles 
situated randomly in a cubic, 3D sample 1.6/^ on a side with particle numbers ranging from 100 
to 937. The material parameters used correspond to ZnO and Ag, over the wavelength range 
300nm < A < 400nm. The particle diameter considered for ZnO is 50nm and for Ag 20nm. The 
eigenvalue spectrum of the scattering matrix is presented and the magnitude of the lowest eigenvalue 
is used to estimate the extent of light localization. The scattering problem was solved for a unit 
amplitude incident plane wave. The energy densities and losses, both total and scattering, were 
calculated. No evidence of light localization was found, neither in the existence of a zero eigenvalue 
nor in any low lying isolated resonance, although the lowest eigenvalues for Ag are an order of 
magnitude smaller than those for ZnO. Energy density calculations show that energy is mostly 
confined to the region near and in the scattering particles, and that little energy concentrates in the 
fields between particles. The energy-energy correlation function shows a peak at the mean particle 
separation in a universal manner. Mixtures of Ag and ZnO particles show isolated spectra with 
little interaction consistent with this energy localization. Gain was added to the ZnO system by 
adding a resonant term having Lorentzian line shape with a negative coefficient, which drives the 
imaginary part of the dielectric constant negative corresponding to gain. With gain in the system, 
the losses can be compensated; however, none of the eigenvalues moves toward zero, indicating that 
localization and loss compensation (which leads to a lasing threshold) are not connected. These 
results indicate that similar calculations with a larger number of scattering sites would still be of 
interest. This, however, will require very different numerical techniques than have hitherto been 
applied. 

PACS numbers: 42.25.Dd;41.20.Jb 



I. INTRODUCTION 



Light localization due to coherent multiple scattering 
and interference in random systems has been the sub- 
ject of much speculation and investigation since the late 
1990's. This speculation arises in connection with ex- 
perimental studies of laser action in random media by 
Lawandy et al. . Their results indicated that an earlier 
diffusion transport theory of the lasing phenomenon by 
Letokhov^] was inadequate to explain the low threshold 
activation energy for laser action. Light localization due 
to interference was suggested by Wiersma et al.Q- Ex- 
perimental evidence was presented, although more recent 
work has called this into question and indicates an impor- 
tant roll for absorption Also, this later experimental 
work, using microwaves, seems to rule out localization 
in 3D. Localization was suggested in an analogy with 
electrons in 3D random systems which has been shown 
will localize|3 . The surmised localization condition for 
light is that the mean free path due to scattering Iscatt 
be approximately equal to \/2tt and that the scattering 
be "strong". This is the loffe-Regel 7] condition for lo- 
calization with a random array of scatterers in 3D. One 
should be able to check this with a computer simulation. 
An unambiguous criterion for a local state to exist is the 
vanishing or near vanishing of any of the eigenvalues of 



the scattering operator. Near localization would manifest 
itself as a sharp separated resonance near zero. Trends in 
the variance of the eigenvalue spectrum with the number 
of scattering sites and frequency could also measure the 
extent and onset of localizationQ- Earlier simulations of 
light in 3D scattering 9] have used no more than 100 scat- 
tering sites. The eigenvalue spectra of these systems did 
not show any appreciable localization. Since the absence 
of localization may be a consequence of system size, we 
are motivated to repeat these attempts with many more 
scattering sites and to examine the scattering operator 
eigenvalue spectrum and also the energy density distri- 
bution when incident light is present; this should give 
some picture of the extent of such localization. 

Accordingly, we present results here for systems with 
up to 900 scattering sites which seems to be near the up- 
per limit for direct numerical methods that still yield con- 
vergent results. In this paper, we use material parameters 
appropriate to ZnO for dielectrics and for Ag for metals 
and confine ourselves to the 300rim < A < 400nm wave- 
length range. These particle sizes are consistent with the 
long wavelength limit. 

Further, there have been reports of induced lasing in 
samples of random ZnO particles in a powder JOJ llj and 
a calculation of a lasing threshold in a model with up to 
900 scattering sites Light localization has been cited 
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as a mechanism that would provide the feedback needed 
for lasing action in these 3D systems. The connection of 
localization and lasing threshold for ID systems seems to 
be well established but it has been questioned for 3D 
systems ^14] 15]. Hence an investigation of the connec- 
tion of this lasing threshold with localization is therefore 
warranted. 

The theoretical formulation of the light scattering 
problem is presented in Sect. II, whose details are rel- 
egated to an Appendix; together these give a review of 
the vector field scattering theory employed here. Its gen- 
eralizations are readily apparent and will be published 
elsewhere. The material model to which this theory is 
applied is presented in Sect. Ill along with the formu- 
lation of quantities for which numerical results are pre- 
sented. Eigenvalue and scattering results are described 
in Sect. IV. In Sect. V nonradiative gain is considered. 
The conclusions which can be gleaned from this work are 
stated in Sect. VI. 



II. THEORETICAL FORMULATION 



scattering method using a free particle propagator or 
Green's function appropriate to photons and scattering 
centers, randomly situated, that represent the medium. 
This formulation is derived in a very general fashion in 
the Appendix. 

The solution for the electric field has the scattering 
equation form: 



E{f)^Eo{r)+ I d^r' G {f~f') ■ k^inPif") (1) 



where r is the source location, r the observation position 
and we define a tensor Green's function as 



G{r-r')= Ij +^\g(\r-f'\). (2) 



The clearest formulation of the light scattering prob- 
lem in random media is probably given by a multiple 



Expanded, Eq.(|2l becomes 
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is the scalar retarded Green's function. In the usual con- 
tinuum formulation of the light scattering problem |l6j| . 
the polarizability is identified as 



and 



AttP (r) = Ae (f) E (r) 



Ae(rl = (e(rO-eo) 
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Here e (r) is the position dependent dielectric constant 
and Co = 1 in Gaussian units. We assume here linear 
response. Other, higher order, induced moments may 
not respond as simply, but we will stick to dipoles as the 
elementary sources. It should be pointed out here that 
in the continuum formulation, where a spatial varying 
dielectric function is used, the dielectric function could 
have a general tensor form. Other boundary conditions 
can be easily applied. For example, there are many 
formulations for implementing periodic boundary condi- 
tions for the scalar Green's function [l^ and carrying this 



over to the tensor case is quite easy as the tensor Green's 
function is obtained by local operations, i.e. derivatives 
of the scalar Green's function, g (R) . 

Here we consider a random medium consisting of N dis- 
crete grains. We shall take the typical dimension of all 
grains to be <g; A, where A is the light wavelength and 
i : 1 ^ N. The integral in Eq.(^ now becomes a sum of 
integrals, one over each grain. We use fi and rj to label 
the location of the centers of the discrete grains i and j, 
respectively. For all integrals over grains where ^ rj, 

Eq. (O can be used and the E (fj ) taken as constant over 
the grain as is appropriate in the long wavelength limit. 
For the case where — fj , the integral over this "diago- 
nal grain" is singular and the integral must be considered 
as a principal value integral. The singularity is extracted 
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where L is due to an exclusion volume at the center of 
the diagonal grain. This singularity is due to taking 
both V's in the definition of Eq.Q under the integral 
sign in the vector scattering equation, Eq.(^. Mathe- 
matically this exclusion or depolarization volume shape 
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FIG. 1: The model sample 1.6/i on a side showing 900 randomly distributed 50nm diameter ZnO spheres. In this case the grid 
on which they are distributed is cubic with 50nm spacing. 



depends on the shape of the grain and must be needle 
shaped if the grain is so shaped etc. This shape depen- 
dence follows from a choice in the order of taking lim- 
its in defining the principle value. These considerations 
are warranted due to the scale invariance of Laplace's 
equation. This singularity has been extensively studied 
by many authors |18j.,.19,j and its extraction is central to 
the derivation of the Lorentz-Lorenz equation when going 
to the continuum limit inside a scattering medium [20| . 



Such mathematical considerations seem non-physical and 
the proper limiting proceedure may be application de- 
pendent. Lorentz originally assumed spherical or cubic 
symmetry plj]. although he indicated the possible exis- 
tence of other terms. In what follows, assume a spherical 

diagonal grain; with spherical symmetry, L— -j- Ex- 
tracting the singularity and integrating over the diagonal 
grain, Eq.QJ now takes the form 
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where Vg^ is the volume of the grain at site j, and Di, 
the diagonal term for site i, given by 



A = i-Ae in) 
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where is now the radius of the spherical diagonal grain. 
This result has been obtained before however, its 
significance is further elucidated by expanding to third 
order in kai. Keeping terms to that order and keeping 
only the first term in the rhs of Eq.©, we have for the 
field inside of a grain: 
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where we have used Eq.© and dropped the site index. 
Eauation (|10|) is now recognizable as the Lorentz correc- 
tion factor to the incident field for the field inside a grain, 
but, with additional dynamic terms. The terms propor- 
tional to fc^ and k^ are called the dynamic depolarization 
and radiation-damping correction , respectively. 

If the diagonal grain is left out and we set Di = 1, this 
Green's function method becomes identical with what is 
known as the coupled dipole method in which vari- 
ous larger grain geometries are represented by clumping 
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the individual dipoles. Then, a very large matrix lin- 
ear scattering problem must be solved with its attendant 
convergence problems, especially if metal grains are mod- 
eled, in which case, one must account for plasmons and 
other dynamic effects due to source kinetic energy. The 
coupled dipole method has been used to check the local 
field approximation 
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used in deriving the Lorentz-Lorenz formula indi- 
cating that, on the average, it is numerically feasible to 
recover the continuum results from the coupled dipole 
method. 

Including the diagonal grain allows us to extend the 
long wavelength approximation to larger particle sizes 
|24l | and to various shapes, although here we restrict our- 
selves to spheres. The scattering equation, Eq.©, can 
be written as 



N 



Y,V{n,rj)-E{r,)^Eo{n), 



(12) 



where 
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For all scatterers of identical species, the matrix V is 
symmetric in all of its indices or, if species vary, it can be 
made symmetric by multiplying both sides by Ae (fi) Vg. . 
However, because the retarded potentials were used, it is 
not Hermitian. Inversion of the complex matrix V, solves 
the scattering problem. Diagonalization of V yields com- 
plex eigenvalues. The scattered field E is decoupled from 
the external field Eg in any mode for which the magni- 
tude of the eigenvalue 0. Such a mode can then be 
considered as rigorously localized. The size of the low- 
est magnitude eigenvalue can be used to characterize the 
extent of localization. 



Power losses were calculated using the optical theorem. 
These are integrals over the cubic sample surface. For the 



power absorbed we have 



Pabs = --^Re 
An 

and the power scattered 



Pscatt = +—Re 
4tt 



Surface 
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where Ei = E — Eo and Hi = H — Ho are the fields with 
the incident light fields subtracted out and Re means the 
real part. The total losses are thus 

Ptotal = Pabs + Pscatt- (16) 

Using Gauss' theorem and Maxwell's equations, these can 
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be transformed into 
ck 

Ptotai = 



and 



iAe*E* ■ Eodv (17) 
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o'"' J vol.scatt. 

These are integrals over the scattering volume, which in 
this case is a sum over the volumes of the N scatter- 
ing grains. Ptotai is of the same form as the extinction 
coefficent |23| ■ Associated with each of the power losses 
there are path lengths, labs,lscatt and Itotai which are cal- 
culated generically as 
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This is essentially the reciprocal of, the crossection as- 
sociated with P times the scatterer density, where we 
assume one species and one grain size for simplicity. For 
example, if we neglect multiple scattering and use the 
familiar 
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for unit incident field amplitude, we can obtain a length 
if'catt'^ ■ This power flow and length leave out multiple 
scattering effects. 



Field energy densities were calculated using 

E{t^ ■ D*{r) 



WE{r) = Re- 
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for the energy density in the electric field, and for the 
magnetic field energy density 



WH{r) = Re 



H{f) ■ H*{r) 



(23) 



where D = eE, for grid sites occupied by dielectric scat- 
tering grains, Re e > 1. For metals, this is not the correct 
formulation of the grain energy as Re e < and so this 
contribution is left out which should make no difference 
for the comparison purposes needed here. The total en- 
ergy contained within the cubic sample is obtained by 
integrating these densities over the sample cube 



Wtotal = 
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where w (r) — we (r) + wh (r)- This total does not con- 
tain the far fields and can only be used as a relative mea- 
sure of energy concentration. We also calculate a figure 
of merit in analogy with resonant circuits 
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In order to characterize the energy density distribu- 
tion, we calculated the energy density correlation func- 
tion, defined as 
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The integrals were evaluated as sums over the sample 
lattice points, the result smoothed and the R deriva- 
tive taken by using a Savitzky-Golay 4 point filter'2^. 
In this calculation, the particle dimensions are equal to 
or smaller then the grid spacing so the contributions to 
the integrals will ignore correlation in the same particle. 
The constant obtained from one particle will be obscured 
when the results normalized by their maximum value. 
Also, the contribution to the sums is limited by the ex- 
tent of the lattice grid; nevertheless, Eq.(j2Hll gives some 
measure of the extent and position of energy localization. 



III. MODEL AND CALCULATIONAL DETAILS 

The model used in generating most of our results con- 
sists of a cubic sample 1.6 microns on a side with grid 
points inside this sample. Using a grid, though not 
strictly necessary, enables the fields to be evaluated out- 
side of the grains more easily. Having 32 grid divisions 
per side produces 32^ — 32768 sites in the sample. This 
means that the grid interval is 50 nm. On a random 
set of grid points, we place our small spherical scatter- 
ing grains with diameters consistent with the long wave- 
length limit. We shall examine grains of ZnO with 50nm 
grain diameters and grains of Ag with 20nm grain diam- 
eters. Figure Ogives an example of the sample with N = 
900, 50nm diameter ZnO spherical grains. When consid- 
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FIG. 2: (a) The dielectric constants of ZnO(reaI part: long- 
dashed, imaginary part: dot-dashed) and Ag(real part: solid, 
imaginary part: dashed) vs A, for SOOrzm < A < 400nm. 
(b) The classical penetration depth for ZnO (long-dashed) and 
Ag(solid) in nm vs A, for 300nm < A < 400nm. 
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FIG. 3: (a) The magnitude of the average eigenvalue and 
the standard deviation (cr) of the eigenvalue magnitude dis- 
tribution vs A for 50nm diameter ZnO scattering grains. The 
average is approximately the same for any number of scatter- 
ers. (b) The standard deviations are given for the number of 
scattering grains in the sample: N = 100, 300, 500, 700 and 
900 (see legend). 



ering light scattering, take the incident beam of light to 
be a plane wave of fixed (unit) amplitude and phase going 
in the z direction with E along x. The range of wave- 
length considered is SOOnm < A < 400nm. Note that the 
plasma resonance of bulk Ag is at about 324 nm and that 
the excitonic emission wavelength that leads to lasing in 
semiconducting wurtzite structure ZnO corresponds to 
385nm, which is also close to the band gap for this ma- 
terial. The ZnO scattering grains have a homogeneous 
dielectric constant whose frequency dependence is shown 
in Fig. [21a) 26J. Both the real and imaginary parts are 
given. Also shown is similar data for Ag [23. Figure|21[b) 
gives the calculated "penetration depth" corresponding 
to these dielectric constants 28J . 

To find the lowest eigenvalues of the scattering matrix 
"P, the scattering problem was solved using the conju- 
gate gradient algorithm _29J with the fiduciary incident 
field described above. No transforms were done as there 
is no benefit to be gained from either the fast Fourier or 
wavelet transforms when the scattering centers are ran- 
domly distributed and there are no discernable internal 
boundaries. This may not be true, however, for a very 



large number of scattering centers, N ^ 900, whereas 
N = 900 is as large a number as could be treated reli- 
ably in the present calculations. 



IV. EIGENVALUE AND SCATTERING 
RESULTS 

For the model described above, using 50nm diameter 
ZnO spherical scattering grains, the magnitudes of the 
eigenvalues of the scattering operator given by Ea. l|13() 
were calculated. The results are summarized in Fig.(|21l. 
The average eigenvalue magnitude is the same for any 
number of scattering grains but the spread estimated by 
the standard deviation, cr, increases with scattering num- 
ber and thus the tail of the distribution extends more to- 
ward zero. There is a slight shift of the average to higher 
values as the wavelength increases. 

The scattering equation, Eq.(|HIl, was solved using an 
incident plane wave of unit amplitude impinging on the 
sample volume from the z direction (see Fig.(^l). With 
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FIG. 4: The losses for 50nm diameter ZnO scattering grains 
for a unit amplitude incident plane wave impinging on the 
model sample Fig.( ^ from the negative z direction vs A: 
(a) gives the total losses Ptotai (black heavy curves) and the 
absorption losses Pabs (gray lighter curves), the number of 
scattering grains in the sample: N = 100, 300, 500, 700 and 
900 (see legend) . (b) gives the scattering losses only Pscatt for 
the corresponding number of scatterers as in (a) above. 
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FIG. 5: The figure of merit Q defined in analogy with that 
calculated for resonant circuits for 50 nm diameter ZnO scat- 
tering grains vs A and for the same range of scattering grain 
numbers considered above. 
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FIG. 6: The energy-energy correlation function C, for ZnO 
A = 385nm and A'^ = 900. The y axis is normalized to the 
maximum correlation and the x axis is normalized to the cal- 
culated average distance between scattering centers lave- This 
curve appears to be universal for any number of scattering 
grains and at any wavelength; i.e. the curves for any wave- 
length, scattering strength and particle number are practi- 
cally identical and completely overlap. The correlation func- 
tion C, calculated for Ag also yields the same curve and over- 
laps the one given above. 



the solution of the scattering problem, the losses and 
field energies were calculated and the results for ZnO 
are summarized in Fig.(0J. Part (a) of this figure gives 
the total and absorptive losses and part (b) gives the 
scattering losses only. As expected, there is a decrease 
in these losses as the wavelength increases. Also the 
losses due to scattering increase significantly as the par- 
ticle number goes up. In ZnO, the process happens in a 
very monotonic manner. These losses can also be ana- 
lyzed in terms of lengths using Ea. H19() . For example, for 
N = 900, kotai « O.lOA at A = 300(nm) ^ kotai « 0.20A 
at A = 400(nm) and labs ~ 0.22A 0.55A, whereas 
hcatt ~ 0.15A 0.21A. The scattering length to wave- 
length ratios increase monotonically with incident wave- 
length as expected. Note that, even for the largest num- 
ber of scattering sites that can be used in this direct cal- 
culation, Iscatt is still slightly larger than the loffe-Regel 
value of « .16A over most of the wavelength range con- 
sidered, and this is with the average distance between 
scattering sites of lav ~ 2.9A for N = 900. It seems, 
that, for realistic materials parameters, the loffe-Regel 
limit is almost obtained in a direct calculation, however, 
no eigenvalues approach zero. 

The relative figure of merit, or Q, for ZnO is given in 
Fig-O, and shows that Q decreases with an increase of 
scattering number. This means that losses dominate and 
that there is no tendency to store energy. The energy- 
energy correlation function, which is defined in Ea. (|26|l . 
is given for ZnO in Fig.(|HJl. The solid lines are for a 
wavelength of 385nm and N = 900. These results are 
normalized to the maximum correlation value on the y 
axis and to the mean distance between scattering centers 
lave on the X axis. There is a broad 
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FIG. 7: (a) The magnitude of the average eigenvalue and 
the varience (cr) of the eigenvalue magnitude distribution vs 
A for 20nm diameter Ag scattering grains. The average is the 
same for any number scatterers. (b) The standard deviations 
are given by the dashed curves for the different number of 
scattering grains in the sample: N = 100, 300, 500, 700 and 
900 (see legend). 

as expected if the energy is localized in and around the 
grains. Note that, as explained above, there is no same 
particle correlation due to the small particle size in rela- 
tion to the grid spacing. There are no peaks off of this 
maximum and the curve appears to be universal, i.e., it is 
practically the same for any number of scattering grains 
and at any wavelength and for any scattering strength. 
This means that the curves for any values of these pa- 
rameters completely overlap and are identical to Fig.(|51). 
It appears that this kind of energy- energy correlation is 
quite general for point particle scattering. 

Results analogous to those for ZnO are given for 20nm 
diameter spherical grains of Ag in Figs.( [T]), ( ISJ and 
(El). Fig. ( [71) for Ag eigenvalues is to be compared to 
Fig. (13) for ZnO eigenvalues. These results can be cor- 
related with the dielectric properties shown in Fig.(|2I). 
For, Ag the scattering strength is resonantly large when 
the real part of the dielectric constant is near -2 which is 
the Lorentz resonance and this yields the lowest eigenval- 
ues and a dip in the standard deviation. The sharpest dip 
in the eigenvalues, however, occurs where the imaginary 
part of the dielectric constant is smallest. These effects 
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FIG. 8: The losses for 20nm diameter Ag scattering grains 
for a unit amplitude incident plane wave impinging on the 
model sample Fig.(0 from the negative z direction vs A: (a) 
the total losses Ptotai (black heavy curves) and the absorption 
losses Pabs (gray lighter curves), the curves corresponding to 
the number of scattering grains in the sample: N = 100, 300, 
500, 700 and 900 (see legend), (b) the scattering losses only 
Pscatt for the number of scatterers as in (a) above. 
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FIG. 9: The figure of merit Q defined in analogy with that 
calculated for resonant circuits for 20 nm diameter Ag scat- 
tering grains vs A and for the same range of scattering grain 
number considered above. 
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become more pronounced as the number of scatterers in- 
creases. By contrast the ZnO results show no such struc- 
ture and only monotonic trends. The eigenvalues for Ag 
are lower than those for ZnO and the spread is much nar- 
rower. These effects are dramatically illustrated in the 
scattering losses. Fig.(|Hll shows the results for Ag and 
compares to Fig.(^ for ZnO. Again the results for ZnO 
show nothing remarkable whereas the Ag results have a 
large resonant peak and fall to a shallow minimum where 
the imaginary part of the dielectric constant in minimal. 
The length results mirror the structure of the power loss 
results in an inverse fashion, however, Iscatt/^ on a scale 
one hundred times larger then either hotai/^ or lats/^- 
Even at the Lorentz resonance Iscatt ~ 1.0 A, however, 
at A = 324(nm), where the imaginary part of e is least, 
Iscatt ~ 600A. It is evident that length reduction due 
to losses is significant and that absorbtion accounts for 
most these losses. 

The Q results for Ag are given in Fig.(|5I) and are again 
dominated by the scattering losses, though it is larger 
than that of ZnO, shown in Fig.([SJ), there doesn't appear 
to be any tendency to store energy since it also decreases 
with increasing scattering number. The energy-energy 
correlation corresponds to the universal curve given in 
Fig.( EJ. This means that all the calculations for Ag 
with various grain numbers and wavelengths result in the 
same correlation curve C, which completely overlaps the 
one given in Fig.(|Sl). 

The eigenvalue magnitude spectra for ZnO and Ag are 
widely separated, as shown in Fig.f llUfl fa'l which gives 
results for 300, 50nm diameter ZnO scattering grains 
and 20nm diameter Ag scattering grains at wavelength 
A = 38Gnm graphed on the same scale. The results run in 
value from zero to 2.75 with a bin size of A = 0.05. and 
the graph maximum for the number of values in any bin 
is 200. The Ag results lie in a much narrower range and 
are lower in value than for ZnO. When mixtures are con- 
sidered Fig.f I10|) fb') is obtained. These results are due 
to a mixture of 150 ZnO grains and 150 Ag grains for 
the same wavelength as used in Fig. ( I10|l (a) . One might 
expect that if there were any phase interference between 
light scattering from these different species having widely 
differing scattering strengths, that there would be some 
overlap of the eigenvalue spectra, but as seen in the fig- 
ure, there is none. 



V. EFFECT OF ADDING GAIN 

In order to explore the connection between light local- 
ization and lasing threshold, gain was added to the ZnO 
system by including a resonant term to the background 
dielectric constant: 

Here es is the ZnO backround, kr is the resonant 
wavenumber taken to be 1.632 * 10^(l/m). This corre- 
sponds to A = 385nm which is at the bandgap in ZnO. 



Eigenvalue 
(a) 



IB 



Eigenvalue 
(b) 



FIG. 10: The eigenvalue spectrum of the scattering operator 
given as a histogram of the eigenvalue magnitudes for A = 
380nm and for: (a) 300 scattering grains of 50nm diameter 
ZnO in black and 300 scattering grains of 20nm diameter Ag 
in gray (b) a mixture of 150 grains of ZnO and 150 grains of 
Ag in dark gray. The eigenvalue range is from to 2.75 with 
a bin size of A = 0.05. The y axis is truncated at a value of 
200 for both (a) and (b). 
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wavelength A (nm) 

FIG. 11: The dielectric constant vs A for ZnO where an 
additional resonant term has been added to the backround 
constant for ZnO shown in Fig.J^J. The real part is in black 
and the imaginary part is in gray. The resonance is located at 
A — 385nm with a width of 7 = 5 * 10^m~^ corresponding to 
the ZnO gap and pre- lasing lineshape given by Cao|lll|. Gain 
levels of 0.0, -1.0, -2.0, -3.0, -4.0 x lO" are shown by the 
dashed lines. The imaginary part of the dielectric constant 
becomes more negative as the gain increases. 
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FIG. 12: The losses for 500 ZnO scattering grains, 50rim in 
diameter and a unit amplitude incident plane wave impinging 
on the sample from the negative z direction vs A: (a) the 
total losses Ptotai (black curves) and the absorption losses 
Pabs (gray curves) and (b) the scattering losses only Pscatt- 
The families of loss curves give results for gains of -2.0, -3.0, - 
4.0 X 10" with the curves dipping progressively more negative 
with increasing gain magnitude. 
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wavelength A ( nm) 
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-2.0*10-^13 
-3.0*10^13 
-4.0*10^13 



FIG. 13: (a) The magnitude of the average eigenvalue of 
the eigenvalue magnitude distribution vs A for ZnO scatter- 
ing sample given above for the corresponding gain levels, (b) 
The standard deviations [a) are given. We note that the distr- 
bution of eigenvalue magnitudes changes somewhat; however, 
the lowest eigenvalue magnitude does not move toward zero. 



The width 7 = 5* 10^(l/m) is taken from the pre-lasing 
lineshape given by Cao et al.[ll|- The gain is added 
by making a negative, simulating a population inver- 
sion. Gain levels a = 0.0,-1.0,-2.0,-3.0,-4.0 aU x 
10^^(l/m)^ were used. The dielectric constant for these 
values in the range of wavelengths 340nm — > 400nm is 
shown in Fig.( lll|l . As gain is added the real part shows 
enhanced scattering on the higher energy side of the res- 
onance and suppressed scattering on the lower energy 
side. The imaginary part below the x-axis and becomes 
more negative as more gain is added. To show that losses 
can be compensated for, scattering results for these gains 
are given in Fig. ( I12|l (a) and (b) . These indicate the 
losses for an incident planewave impinging on the sam- 
ple from the negative z direction, as previously discussed, 
for our model. In this case, however, results for 500 ZnO 
scattering grains are given but now with the indicated 
gains included in the dielectric constant. Note that total 
losses can be driven negative. Scattering losses are also 
affected, but of course remain positive. The only length. 



that is non-singular and makes sense for this system, is 
Iscatt which, mirror the Pscatt results in an inverse fash- 
ion. For the largest gain, Iscatt ~ 0.55A at A = 383(nm), 
its lowest value, below the resonance and Iscatt ~ 2. OA 
at A = 394(rtm), its highest value, above it. Results are 
shown for gains -2.0, -3.0, -4.0a:10"(l/m)2. Finally, 
Fig. ( [T^ gives (a) the average eigenvalue magnitude and 
(b) the standard deviations in the relevant frequency 
range for these gain levels. The spectrum is affected by 
the gain, the average dipping lower on the low energy 
side of the resonance, even though the scattering would 
seem to be somewhat suppressed there as indicated by 
the dielectric constant results. In any case, none of the 
eigenvalue magnitudes were driven toward zero. This 
indicates that, at least in 3D systems, localization and 
lasing threshold are not connected in any direct manner 
and that the interpretation of random laser experiments 
should be re-examined. 
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VI. CONCLUSIONS 

It can be seen from the above results that there is very 
httle energy concentrated between scattering particles, 
but that most of the energy is localized around and in 
the grains. Energy concentration due to phase interfer- 
ence does not seem to happen for light in a 3D random 
scattering medium, at least for the system dimensions, 
number of scattering centers and material parameters 
considered here. There is no zero eigenvalue represent- 
ing a local state and no isolated resonance near zero that 
appears for a critical number of scattering particles. The 
loffe-Regel condition is not quite met with realistic ma- 
terial parameters in this direct calculation. However, the 
tendency to localization with randomness or with an in- 
crease in the number and/or strength of the scatterers 
seeems to be gradual and continuous with no indication 
of any sharp transition. The energy-energy correlation 
calculations, which show a broad peak at the mean par- 
ticle separation, and the fact that the curve seems to be 
universal, buttress this conclusion. Also, there is little 
interaction between the spectra of scatterers of different 
strengths, as is indicated in the spectra of Ag and ZnO 
mixtures. 

Adding gain to such a system of random scatterers, 
even enough to completely compensate for losses, does 
not cause the magnitude of any of the eigenstates to move 
to zero. This is in spite of the loffe-Regel condition not 
quite being met. The lasing threshold can be consid- 
ered as the point where the gain just compensates the 
losses. Thus it would seem that localization and lasing 
threshold have little connection. This means that the in- 
terpretation of the mechanisms of random lasing should 
be revisited. While we have used an order of magnitude 
more scattering sites that was done previously, it is clear 
that a still larger number of scattering sites would be of 
interest in order to check for localization in realistic 3D 
random system. Experimental occupied volume fractions 
of 50% are auotedfllj and, in the terms used here, that 
amounts to about N « 30000. For this simulation, how- 
ever, very different numerical techniques will have to be 
used than have hitherto been applied. 
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APPENDIX A: 



DERIVATION OF GREEN'S 
FUNCTION 



Using the retarded potentials: 

1 f j{r' ,t')(Pr' 



and 



A{r,t) 



R 



p[r' ,t')(Pr' 



R 



(Al) 



(A2) 



where p and j are the source charge and current densities, 
R = \f — f I, and t = (t — -^) is the retarded time (i.e. 

the time that it takes light to travel from r to r^, the 
Lorentz gauge is chosen so that 



V- =0. 

c ot 



(A3) 



The continuity equation for the source terms must be 
satisfied, i.e. 



^' ^ I dp „ 



(A4) 



and this is manifestly true if the source terms are given 
by 



dP 



dQ 



J = ^+cV xM-V •^ + ... (A5) 



dt' 



and 



-V • F + V V 



(A6) 



which is also consistent with the Lorentz condition and 
the representation of the potentials given in Eqs. pri 
IA2|) . Here P is the electric polarization density, M the 

magnetic polarization density, Q the electric quadrupole 
density, etc. The derivatives here are taken with respect 
to the source location r and time t' . The elementary 
charges can be grouped at a point into moments of any 
order to represent the material medium both statically 
and dynamically; however, here we shall consider only 
electric dipole terms for simplicity and to illustrate the 
technique. 

We now make explicit the derivatives with respect to 
the observer position contained in the retarded time and 
write 
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V •[P] = [V -P] 



dP 
'dt 



dP 
'dt 



R_ 



(A7) 



r 



where R = r — r and [ ] denotes evaluation at retarded time, thus 
I 



V P 
R 









• 




v' • 






P 


■ R 


p 















R 



cR'^ 



-cfr'. 



(A8) 



Also note that 
\P 



P 



R 



R 



P 



the Green's function formulation of the scattering prob- 
lem and with the coupled dipole approach 24J we 
shall assume a r 
(A9) this assumption: 



so that 



V • 



P 



P 



R V 



P 



R 



i?3 



R 



(AlO) 



Therefore, we finally have: 



P 



R 



i?3 



-(fr 



• 

p 


■R 


cl 


'?2 



(All) 



where P indicates a derivative with respect to time. In 
a similar way we obtain 



A(f, t) 



P 



c . R 



The fields are given by: 



and 



I OA 
c ot 



H = \7 X A. 



(A12) 



(A13) 



(A14) 



We need only solve the equation for the E field because 
from it we can obtain H. 

Up to this point, we have made no assumption about 
the time dependence but in order to make contact with 



P (t) = -iiuP (w) e" 



Pit) 



(A15) 



(A16) 



Letting k — — we have 



P{t) 



= -ickPe'^''^-'^'l 



(A17) 



Note that the important e**^^ phase term is due to evalu- 
ation at the retarded time. The potentials are now given 



as 



A{f, t) = -ik / d'^r 



P r 



J{kR-u}t) 



R 



so 



IdA 
'c~dt 



-"-11 = e I d^' 



,P{r') 



,i{kR-Ljt) 



R 



(A18) 



(A19) 



and 



(A20) 

Using Ea. (|A12p and Ea. HAll|) . Factoring out the e"*"* 
time dependence, we have for the E field: 
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E{r)^Eo{r) + 



,P{r') 



ikR 



R 



(fir 



ik 

'r 



1 



RP 



(0 



■ikR 



R 



(A21) 



where we have added an external "incident" field which 
comes from other sources not given by the potentials con- 
sidered above; this is the inhomogeneous term in the light 
scattering equation. The H field can be obtained from 



H{r) = V X E 

k 



(A22) 



which is Faraday's law using our monochromatic assump- 
tion. Note that the V term drops out when the curl 
in the above equation is applied to the E field given in 
Ea. HA21|l . The scalar retarded Green's function is de- 



fined as 



so that 



9{R) 



JkR 



AttR 



Vg{R)^g{R) ['-^^±]R. 



(A23) 



(A24) 



Using this relation Eq. (|A21|I becomes Eq. (QJ in the main 
text above. 
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